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Abstract 

In this paper we describe the implementation of the QED process 77 77 through a fermion loop 
into the framework of SANC system. The computations of this process takes into account non-zero 
mass of loop-fermion. We briefly describe additional precomputation modules used for calculation 
of massive fermion-box diagrams. We present the covariant and helicity amplitudes for this process 
and also some particular cases of Do and Co Passarino-Veltman functions. Whenever possible, we 
compare the results with those existing in the literature. 
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1 Introduction 



The SANC is a computer system for semi-automatic calculations of realistic and pseudo-observables for 
various processes of elementary particle interactions " from the SM Lagrangian to event distributions" at 
the one-loop precision level for the present and future colliders — TEVATRON, LHC, electron Linear 
Colliders (ISCLC, CLIC), muon factories and others. To learn more about available in SANC processes 
see the description in [1] and look at our home pages at JINR and CERN [2J. 

Light-by-light scattering is one of the most fundamental processes in QED. It proceeds via one-loop 
box diagrams containing charged particles. The first results for the low energy limit of this process 
were obtained by Euler [3J. Then Karplus and Neumann [1] found a solution for QED in general but 
complicated way. The QED cross sections in the high energy limit were calculated by Ahiezer [5]. 
Nowadays there are computations for 77 — > 77 process in the electorweak Standard Model by Bohm [5] 
and even for two-loop QCD and QED corrections by Bern [7] . 

In this paper we describe the implementation of the QED process 77 — > 77 through fermion loop 
and corresponding precomputation block into the framework of SANC system. The computations of this 
process take into account non-zero mass of loop-fermions. 

The paper is organized as follows: 

First we discuss some notations and common expression for cross section in section [2. II 
In section [2.21 we discuss diagrams for 77 — > 77 process and covariant amplitude tensor structure. 
In section 12.31 one obtains compact form for this structure. The idea of form factors is described in 
section 12.41 

The helicities amplitudes approach and their expressions for light-by-light scattering in general (mas- 
sive) and in limiting (massless) cases are listed in section [2751 

In section |3] we shortly describe precomputation strategy [Tj and the place of this process on the SANC 
tree. 

At last in section 3] one can find the result in the limiting case and comparisons with those existing 
in the literature. 

Additionally, in Appendix section [5] we list answers for particular cases of Dq, Cq and Bq Passarino- 
Veltman (PV) functions [8] (see also the book P). Finally, we present a table of integrals over the 
scattering angle, which are needed for calculation of light-by-light scattering through massive and massless 
loop fermions. 

2 Light-by-light scattering process 
2.1 Notation, cross section 

The 4-momenta of incoming photons are denoted by ki and /c2, of the outgoing ones — by and k^. 
The amplitudes are calculated for scattering of real photons, that is fc^ = , A;| = , A;| = , fc| = 0. 
The 4-momentum conservation law reads: 



fci H- fc2 — fcs — fc4 = . 



(1) 



The 



Mandelstam variables are0 



u ~ — 



s = — 



(fci + k2 f = -2ki ■ k2 

(fci - fc4)2 = 2ki ■ ki , 



'2 , 



i--(fci-A:3)2 
s+t+u^Q 



2fci • k 



(2) 



For the 2 2 77 ^ 77 process, the cross section has the form: 



da- 



77-^77 — 




|^dci>(2), 



(3) 



Note, that in SANC we use Pauli metric. 
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where j = 4-^/ (fcifc2)^ , is the flux, A-y^^-y-y is the covariant amphtude (CA) of the process, and d<i>^^) is 
the two body phase space: 



(2) 



(27r)^5 (fci +k2-k3- ki) 



d-^ksS (fc|) d'^k4S (fc|) 
(27r)3 (27r)3 



For the differential cross section one gets: 



da^ 



1 



' d cos 9 , 



'77-77 - ^28^u;2 I^T7-77I 

where ui is the photons energy and 9 — the scattering angle in the center of mass system (CMS). 



(4) 



(5) 



2.2 Covariant amplitude 

The covariant one-loop amplitude corresponds to a result of the straightforward standard calculation of 
all diagrams contributing to a given process at the tree (Born) and one-loop (1-loop) levels. 

The CA is being represented in a certain basis, 
made of strings of Dirac matrices and/or external 
momenta (structures), contracted with polariza- 
tion vectors of vector bosons, e(fc), if any. The 
amplitude also contains kinematical factors and 
coupling constants and is parameterized by a cer- 
tain number of Form Factors (FFs), which are de- 
noted by J-i, in general with an index i labeling 
the corresponding structure. The number of FFs 
is equal to the number of independent structures. 

The 77 77 process in QED appears due 
to non-linear effects of interaction with vacuum, 
so this process has no Born or tree level. Corre- 
sponding diagrams start from the one-loop level 



q + k 




q + k4 



Figure 1: st-channel diagram for 77 — > 77 process 



and in QED there are box diagrams with four internal fermions of equal mass. The number of not identi- 
cal diagrams (or topologies) is equal to six. But three of them differ from another only by the orientation 
of the internal fermionic loop, giving the same contribution or a factor 2 to the amplitude. So, only three 
topologies (st, su and ut channels) remain which are related by simple permutations of external photons 
in the diagram shown in Figure [T] st-channel is shown, su-channel is obtained by k^ <-> k^ and ut-channel 
— by A:2 ^ ^3- 

The full CA of given process for off-shell photons (fc^ei ^ 0) can be written as: 



43 



A 



ot/Sfiv 



(6) 



where e is the electron charge, Q/ is the fraction of charge of loop fermion in units of electron charge. 



Tj '^'^ are tensors defined with an aid of auxiliary strings Tj from the following subsection and J^i are 
FFs, depended on invariants s ,t ,u and also on fermion mass and Passarino-Veltman functions. The 
off-shell process contains 43 basis elements. 



2.3 Strings and basis 

To obtain a compact form of structures of the amplitude we choose 14 auxiliary tensorial strings: 

rf" =k2i,k3p + \t5pi, , 
r^^^k4i3k2f^ + ^uSi 



^1 


= klfjk2a + 




LIU 

^2 


= ^3^1^41/ + 




OiLL 

^4 


— klf_ik4a -\- 


2^^ctfJ- : 


au 
^5 


— kii/k^a + 




^7 


= /ci^ - tu' 


^k2f, , 


^8 


= kii, — ut'' 


^k2i, , 


•7-" 
^10 




'^ksa , 


<1 


— k4fi , T^2 


= kzu 



st 



^k-ifi , 

14 — kia. 



(7) 
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The complete basis T°'^^'^ can be presented in a compact form with an aid of the auxiUary strings Tj : 

^5 — ' 2 ' 9 ~10 ' 6 — ' 3 ~7 ' 10 ' ^7 —~4~S~9' ^8 — ~5 '7 '9 ' 

rpapfiv _ Bn^ a rpa/3iJ.i^ _ BiJ. u a rpaflfiu _ af3 u i, rpafffiu _ uiy 13 a 

-'g ~~6 ~8~10 7 -^28 ~~6 ~12~14; ^11 — ~1 ~7~12 7 12 — ~2 ~9 ~14) 

rpapliv _ JBu a a rpaPfj.}y _ ,y 3 T^a/3^i/ _ qi/ u £) rpaP^}^ _ Bfi i, a 

^13 —~3 ~7~147 -^14 —~4 ~8~13' ^15 " ~5 ~7~13' 16 — ~6 ~8 ~14) 

rpapnu _ al3 u y rpaPfiu _ uu f3 a rpaPtiy _ By u a rpafffiu _ afj,^ B 

^17 —~1 ~11~8 1 -^18 —~2 ~13~10' ^19 — ~3 ~11~10' ^20 — ~4 ~12~9 ' 

rpaf3ij.v _ au^l^ ^1^ rpaPtiv _ B^i „ a rpaPiiu _ af) u „ rpaP^u _ uu B a (o\ 

^21 —~5 ~11~9 ' -^22 —~6 ~12~10i ^23 ~ ~1 ~ll~12i 24 — ~2 ~13~14: \° ) 

rpaptiy _ (By a rpOiPfi}^ _ afj.^ 13 rpa(3 fiu _ av ^t^- J rpOiPfiy _ u (3 a 

^25 ^11^14) ^26 ~^4 ^12^13' ^27 "^5 ^11^13' 10 "^7 ^8 ^9 ^10' 

rpaptiu _ B a rpaPf^u _ n B „ a rpoPfii^ _ n a ^i^ ^0 rpaf3nu _ u B u a 

^29 —^7^8^13^14' -^30 —^7^9^12^14' ^31 —^7^10^12^13' ^^32 —^8^9^11^14' 

rpapup _ u^a ^13 rpaf3nv _ B a rpaPfiu _ u „ B a rpafifiu _ u i, a ^0 

-'33 —^8^10^13^11' -^34 —^9^10^11^12' -^35 — ^7^8^9^14' -^36 —^7^8^10^13' 

rpapiiu _ a (3 a T^a/^M'' _ _i>_/3_a rpa/Bfiv _ a v 13 a rpafffiu _ u i, a ^0 

^37 — T^7 ^9 ^10^12' -^38 —^8^9^10^11' -^39 —^11^12^13^10' -^40 —^11^12^14^9' 

rj^OLpfiu _ u J3 a V rpafBtiy _ u 13 a rpal3iLV _ u v 13 a 

^41 — ~11~13~14~8 ' -^42 — ~12~13~14~7 ' -^43 — ~11~12~13~14- 

Thus, one gets a minimal number of tensor structures of the CA, which contains a lot of terms. It can 
be written in an explicit form with an aid of scalar FFs. All masses and other parameters dependences are 
included into these FFs, but tensor structure with Lorenz indices are given by basis ([8]). It is important 
to emphasize that each basis elements x"^^'^ , i = 1 43 is four transversal with respect to each external 
photon: 

T^f^'^'^k^ = Tf'^'kp = Tf'^'k^ = Tf^^k, = . (9) 

2.4 Form Factors 

FFs are scalar coefficients in front of basis structures of the CA — projections of CA to complete basis 
expressions T"^^'^ . They are presented as some combinations of scalar Passarino-Veltman functions Aq, 
Bq, Co, Dq 8;. They do not contain UV poles. 

The number of terms in FFs equals to thousands in the case of non-zero mass of the loop fermion, but 
this number reduces greatly for zero loop fermion mass. Full answer for FFs one can find in the computer 
system SANC on servers [2_. For massless loop fermion the FFs are rather compact: 

_ 4 /3i7r iir in n'^ut Stt^ 1 \ 4 / ut 3 \ 

4 / SiTT 4i7r iir t 1 1 \ , 2 f ut 3 4 1 \ ,n 

+ o -TT^ ^-- + - + — H*-o -— + - + — + — Ht 

3 \ zs^ su s'^ su J 3 \ su J 

4 f 3iTT Ain ITT u 1 1 \ 2 f ut 3 4 1X3 
+ 3 \2s^ ~^'¥'1^^^^7t)"'3 \^ '^^^Vt^^J 

_ 4 /iTT ITT inu 2iTT TT^ 27r2 Stt^ 1 \ 4/1 4 3 \ 
' " 3\Vt^^~~~l^~2^~~~^^l^)^i\s^^Jt^2^)''' 
4 /4i7r Ai-K Zin in 1 3 V 2/1 4 3 1 \ ^2 

3 \ st su 2t'^ u^ su y * 3 \ su t'^ u^ ) 
A (insu 3i7r 1 u l\ 2/1 4 u^ u 3\ 2 

+ 3 \ l^~2^'"Jt"t^^i^) ""slv^ + ^ + F + f^+f^j "' 
^ _ A fin in int 2in n^ 2n^ Sn"^ 1 \ 4/1 4 \ n 

^ 3 \st su u^ u^ 2s2 su Au"^ u^ J 3 ys^ su 2u'^ ) * " 

4 /4i7r Ain iin in 1 3 \ 2/1 4 3 1X2 
^3 \~ ^^"2^"t2+^~^j ""S^v^^^^^^i^j n 

A ( int int^ iin 1 t 1 \ , 2/1 4 t^ t 3 



3 V M'^ 2u'^ ' su u^ ' u^ ) ^ ?> \ s"^ ^ su ^ u'^ u^ ' u^ ) * ^ 
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4 /SiTT 2i7r 2i7r txH^ txH 1 \ 4 / t 



^7 



where 



3 V si su 2s4 2s3 s2 / 3 V 

4 /3i7r 7i7r 2i7r t 2 2 \ 2 / 1^ ^ ^ ^ _ 
'3l2?^2^^^"^"^"^j *"3l^^^"2?"2^"^j * 



4 /3i7r 7z7r 2i7r t 12 



S2 2si f2 ^3 g2 gi-J " 3 \^54 g3 2s2 2si i2 ^ « ' 

4 /iTT ZTT 2i7r 7r2t 37r2 1 \ 4 /2t 



3 V si su u2 s3 4s2 su J 3 



4 /2i7rt 3i7r 2 1 2 

3 \ 1*3 2li2 s2 su u^ J ^ 3 \^s3 ' 4s2 



4 / i 3 
k — — 









t 


3 


7 


+ -3 - 
s3 


2s2 


2su 


t 


3 


7 






2si 


2s2 J 






t 

+ — 


3 ^ 




u-^ 


4^2 



'3U2^s2 sti" 3U3^4s2 2^2'' 



u > 



^ _ 4 /2i7r 2i7r in n^t Stt"^ 3tt^ 1 \ 4 / 2t 3 ^ W ; 

^ 3 \ st su u^ s3 4s2 4su su ) 3 \ s3 2s2 2su^ * " 

4 / iirt 2 1 1\ 4/t 3 3 ^ \ ,2 

4 /3i7r 3i7r 2i7r 2 2 3 \ ^ f ^ ^ ^ 3 1 
~3\2^^2^~1^^1^^7t^^) ""3 (v"^^4^^4^^2^"t2. " 

_ 'zTT in in n'^f n^t n^ 2\ ( t 1 ,,, 
" ^';^ + ^ + ]^"^"2^~4;^^^j+M^ + ^+2^''*'" 



«7r 


in 


nV 


^2t 


^ + 


su 


2s4 


2s3 


SU 


in 
+ — - 


t 

s-^ 


-)'• 

su / 


in 
st 


in 


S3 





, in in in t 1 \ f t 1 1 1 , 

-4 -^ + — + ^ + ^ M*-2 -r + ^ + ^ + — + 



2s2 su 1/2 s3 SU/ ys** s3 s2 su 1*2 

. I in in in u 1 \ , „ / 1^ t 1 1 1 \ ,9 
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Zt = ln(--V /„ = lnf--). (11) 



s / v s 

Also there are equations among six FFs (even for massive case): 

1*2 

— ^5 , ^6 — "^-^7 , T^ — T^. (12) 

The other FFs (namely 11-43) are not zero, but they do not contribute, because of corresponding basis 
structures for on-mass-shell photons satisfy Ward identity: kiei{ki) = 0. 



2.5 Helicity amplitudes 

In SANC we use helicity amplitudes approach. 

In the expression for C A as one can see in subsection 12.21 one has tensor structures and a set of scalar 
FFs. To calculate an observable quantity, such as cross section, one needs to make amplitude square, 
calculate products of Dirac spinors and contract Lorenz indices with polarization vector. In the standard 
approach making amplitude square one gets squares for each diagram and their interferences. This leads 
to a lot of terms. 

In the helicity amplitudes approach we also derive tensor structure and FFs. But the next step is 
a projection to helicity basis and as a result one gets a set of non-interfering amplitudes, since all of 
them are characterized by different set of helicity quantum numbers. In this approach we can distinguish 
calculations of Dirac spinors and contraction of Lorenz indices from calculations of FFs. We can do this 
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before making square of amplitude. So, proceeding in this way, we get a profit on calculations time (less 
amount of terms due to zero interference) and also more clear step-by-step control. 
For the process 77 — > 77 one gets: 



A. 



■77^77 

|2 



spins 



spins 5 



(13) 



spins 



Note, the total number of HAs for this process is equal to 16. This corresponds to different combinations 
of external particles spin projections. For 77 77 processes there are 4 photons with two polarizations 
'+' and '— ', so the total number is 2 • 2 • 2 • 2 = 16. Helicity amplitudes are scalar expressions. They 
result from an application of the procedure TRACEHelicity .pre: 





= 71 




= n 


n-+— 


= n 




= n 


H+ — + 


= n 


H++ — 


= n 



+ = W — +- = 



1 ' 



+ 

-+-^ 
— 



s Tr, — suTy + suTc, — s J-ic 



— 2suTct + s^Txz 



(14) 



Therefore, at this stage we observe five independent HAs, while in the case of zero loop fermion mass 
one gets only four independent HAs which are very compact: 



— — Ji — ++ - -1 -f 

— \ — = 7Y — \ — \- = —1 — in 

"Ha h = W — hH — = ^1 ^ 



{l + iiT) 



s 



t - s 



2m 



t 



(1 + i7r) 



t 



k 



1 _ S^\,2 

2 



1 _ su\ 2 

2 t2 I'' 



(15) 



All the others HAs are equal to -1-1. 

The relations among helicity HAs are due to C, P, T-invariance. Moreover, another relation is due to 
crossing symmetry: 



n+ — + (s, t, u) = 71+-+- (s, u, t) , 
but this fact does not mean reducing of the number of independent HAs. 



(16) 



3 Process 77 77 in the SANC tree 

For boxes the SANC idea of precomputation becomes vitally important Ij. Calculation of some boxes 
for some particular processes takes so much time that an external user should refrain from repeating 
precomputation. Furthermore, the richness of boxes requires a classification. Depending on the type of 
external hues (fermion or boson), we distinguish three large classes of boxes: ////, f fbb and bbbb. 
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^ Root 
t [[f SANC 
9 €3 QED 



In this section we briefly describe modules relevant for 77 — > 77. The sum of contributions of fermionic 
loop boxes form a gauge- invariant and UV-finite subset, which is a consequence of Ward Identity Eq. 

The precomputation file AAAA Box (see SANC tree in 
Figure[2]) contains the sequence of procedures for calcula- 
tion of the covariant amplitude. At this step we suppose, 
that all momenta are incoming (denoted by p's) and pho- 
tons are not on-mass-shell. Therefore, these results can 
be used for other processes which need these parts as 
building blocks. 

When we implement the process 77 —^ 77 (see 4A 
branch), we use this building block several times by re- 
placing incoming momenta p's by corresponding kine- 
matical momenta fc's, and calculate FFs by the mod- 
ule AA->AA (FF), then helicity amplitudes by the mod- 
ule AA->AA (HA) and finally — the differential and total 
process cross section by the module AA->AA (XS). 

Also in SANC system one has an opportunity of send- 
ing the analytical results to numerical evolution [1] . 



Precomputation 
^ □ Self 

□ Vertex 
^ Box 

INI Box 
-F^ IIAA Box 
- ^ AA FF Box 
'— ^ AAAA Box 
Processes 

□ A -> I I 

□ 41 

□ 2I2A 
^ 4A 

AA -> AA (FF) 
AA -> AA (HA) 
¥m AA -> AA (XS) 



□ BN 
lH QCD 



4 Results and comparison 



Figure 2: 77 — > 77 process in the QED tree 



Now let us summarize the result. 

The differential cross section of 77 



77 process in QED has a form: 



dcT^^ — — g 2 ^ ^ l^spinsl d COS . 
spins 



(17) 



where uj is the photons frequency, 9 is the scattering angle in CMS and helicity amplitudes are expressed 
in terms of FFs. All dependences on Mandelstam invariants and loop fermion mass, also from Passarino- 
Veltman functions are included into these FFs. All intermediate results for this process can be found in 
the SANC system on the project servers [2]. The final answer for the total cross section in the massless 
limit after substitutions of helicity amplitudes and angular integration has a form: 



X --y 'y — J. --y 'y 



/108 



- 8^(3) + - 24C(5) 



(18) 



This result was compared with 5 and the complete agreement was found. Also the limit of helicity 
amplitudes was compared separately with [7\ and again full agreement was observed. 

One should emphasize that the obtained building block for box diagram in QED (as in 77 — > 77) is 
the first step in the creation of environment for calculations of similar processes in the Standard Model 
(like 77 77, 77 ^ ZZ [TU]) and in QCD (like gg 77, gg ZZ, gg W+W' etc.) 
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5 Appendix 

To obtain the cross section of the 77 — > 77 process in an analytic form we have to compute in an explicit 
form the master integrals, Bq, Co, Do — scalar PV functions [8], |9j — for a particular set of parameters. 
In Dq and Co functions one can see collinear divergences, but the differential cross section is free of mass 
singularities which completely cancel in the sum of all terms. The Aq and Bo functions contain UV 
divergences, which cancel in the sum of box contributions. In the process of computation we face also a 
problem of the "angular edge" divergences, but they are not physical and also cancel completely. 

5.1 Bq function 

The Bo function for 77 — > 77 process reads: 



Bo{Q^;M,M) 



2 



2-ln ( ^ 1 - /31n 



A* 



where 



and 



4M2 



13+1 
13-lJ ' 



(19) 



= A/2 - ie 

In zero limit of fermion mass: M — > 0, one gets: 

Bo (-s; M, M) ^ T + 2 - 
Bo{-t;M,M) = T + 2- 



In I — I - ZTT 



k + In 



Bo i-u; M, M) = - + 2 - 



(20) 



(21) 



5.2 Co function 

The Co function is: 



Co 



'0 (0, 0, m, M, M) = I dx I dy (q^y - Q^xy + M^) 







After calculations: 



where 



Co(0,0,Q2;M, M,M) 









MB 







(22) 



(23) 



a;i,2--(l±/3). 



(24) 
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For M -> 0: 



Co (0, 0, -s; M, M, M) 

Co (0, 0, -t; M, M, M) 



1 

'2s 
1 

"2t 



m I ) + i-K 



Co(0,0,-n;M,M,M) = - — 

zu 



5.3 Do function 

The Do function for st-channel diagram looks like: 

1 :r y 



In 



(25) 



Dq (0, 0, 0, 0, -s, -t; M, M, M,M) = j dx j dy j dz (txy -{s + t)xz + syz -ty + tz + M^) ^ (26) 
After lengthy calculations: 







Do (0, 0, 0, 0, -s, -t; M, M, M, M) : 



+27r2y dx 



1 

+27ri j dx 



1 

+27ri j dx 




a; + ^1 \ ^ ^ a: + Ai 



a;- A3 



x + A^, 



a; - A3 



a; + As 



■/. . x + A2\ ( . . x + Ai\ 



1 

+27ri y dx 




a;- A3 



-A2 + -4:, . 



X - A2 
~hTA, 



a; + As 



,T - A2 

-A2 - A, 



a; - A3 



a: + As 



In (1 + As) -In (As) 



In + In + In (Ai - As 



+ In (-As - Ai) + In (A2 - As) + In (-As - A2) 



+ In (-As + 1)- In (-As) 



In 



+ lnr-^ 



.4M2/ V4M2 
+ In (As - Ai) + In (A2 + A3) + In (As - A2) 



In (Ai + As) 



+Li2 



-Ai - As 



-A1+A3I 



+ Li2 



1 + A3 
A1-A3 



Lio 



1-As \ 
Ai+AsJ 
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10 



+1 



■4 + 8a; - 12x^ + 8x^] dx 



331 

144' 



-1 
+1 



In (a;) In (1 - a;) 4 - 12a; + 12a;^ + tt^ (1 - 2xf + An^x'^ (1 - x 



-1 

+1 



1 / , 77397r2 7tt^ 



In^ (1 -a;) (8- 12a;+ 12x^ + 7:2 (3-4x + 8x2 -8x3 + 4a;^)) da; = - ( 125+ ., 



1 2 



dx + y In (1 - x) 



2 3 



- 4 + 4x + TT^ I 3 + 4x - 6x^ + 4x^ 



yin^a-x 

-1 
+1 

j In (1 - x) 

-1 
+1 

J TT^ (3 - 2x + 2x2) + 8 + tt'' Q - X + 2x2 - 2x3 + x''^ 



+1 

/ 



1 1 



X^ X 



dx=iri-^-^ 



3 r 



dx = l + 



llTT^ 47r* 



12 3 



77r4 87r2 
dx = -— + —r- + i 
60 3 



-1 
+1 



+1 



^4 



In-* (1 - x) ^dx = 24C (5) , j In^ (1 - x) ^dx = ~, 



X 15' 



-1 
+1 



-1 
+1 



/ In^ (1 - x) ^dx = ^, / In^ (1 - x) -dx = 2C (3) . 

7 X o J X 
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